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INTRODUCTION 


In  this  report  aerosol  characterization  measurements  are  analyzed  for  predic¬ 
tion  of  countermeasure  aerosol  effects  on  high-energy  laser  (HEL)  beams. 
Although  the  general  principles  are  applicable  for  various  beam  parameters, 
our  calculations  will  be  oriented  toward  those  of  a  pulsed  CO  laser  (10.6pm 
wavelength)  with  a  pulse  length  of  about  3ps  and  a  total  beam  energy  of  about 
200  d.  The  focus  of  this  effort  will  be  the  effect  due  to  evaporative  clear¬ 
ing.  The  basis  for  thermal  blooming1  (atmospheric  lensing  due  to  thermally 
created  density  gradients  in  the  beam),  the  other  major  nonlinear  effect,  will 
be  established  by  describing  heating  of  the  ambient  gases  and  particles. 
Threshold  calculations  show  that  the  3ps  pulse  of  these  tests  is  short  for 
strong  blooming  effects. 

These  measurements  were  performed  in  conjunction  with  a  test  of  an  HEL  having 
the  parameters  described  above.  The  beam  traversed  a  trench  approximately 
100  m  long.  The  trench  was  covered  with  polyethylene  forming  a  relatively 
closed  environment  in  which  the  smoke  was  dispersed  and  then  measured  as  a 
function  of  settling  time.  Most  of  the  characterization  was  performed  at  the 
midportions  of  the  trench  and  assumed  to  approximate  a  spatial  average  since 
the  smoke  mixing,  diffusing,  and  settling  times  were  long  at  1  to  2  h. 


MEASUREMENT  SYSTEMS 

Table  1  lists  the  instrumentation  used  for  the  characterization. 


TABLE  1.  NP  CHARACTERIZATION  I NSTR' MENTATION 


Measurement 

System 

Manufacturer 

Extinction  absorp¬ 
tion,  coefficients 
at  10.6um 

COz  laser 
spectrophone 

Atmospheric  Sciences 
Laboratory  (ASL) 

Particle  counting 
and  sizing 

Active  cavity,  light 
scattering  particle 
spectrometer,  HeNe 

Particle  Measurement 
systems  (Model 
ASASP-X) 

White  light  volume 
scattering  coeffi¬ 
cients 

Nephelomefer 

Meteorological 
Research  Instruments 
(Model  20S0B) 

Relative  humidity  and 
partial  pressure  of 
water  vapor 

Cew-point  hygrometer, 
thermometer 

Edgerton-Germehausen 
and  Grier  (Model 

m) 

Af .  G.  Gebhardt,  1976,  "High  power  laser  propagation,"  Appl  Opt,  15(6) -.1479 
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The  following  paragraphs  describe  instruments  and  discuss  the  need  for  each 
type  of  information  obtained. 


Of  the  absorbing  gaseous  constituents  at  10.6um,  water  vapor  is  expected  to  be 
by  far  the  strongest.  The  very  high  temperatures  and  dew  points  produced  in 
the  trench  by  the  greenhouse  effects  resulted  in  relatively  high  partial 
pressures  and  thus  high  absorption  coefficients.  The  Edgerton-Germehausen  and 
Grier  (FG&G)  condensation  plate  dew-point  hygrometer  was  used  to  obtain  the 
water  vapor  partial  pressures  from  which  the  water  vapor  absorption  coeffi¬ 
cients  were  computed.  Relative  humidity  values  obtained  from  the  dew  point 
and  temperature  data  were  used  to  compute  fractional  water  content  of  the  very 
hygroscopic  smoke  of  burning  phosphorus.  This  computation,  in  turn,  was  used 
to  determine  the  appropriate  complex  refractive  index  and  other  physical 
properties  of  the  smoke.  Although  site-peculiar  gases  are  possible  contribu¬ 
tors,  contributions  to  the  net  10.6pm  absorption  coefficient  due  to  other 
ambient  gases  can  be  reasonably  estimated. 


The  distribution  of  particles  with  respect  to  size  is  important  since  the 
calculations  of  the  nonlinear  effects  are  also  generally  size  dependent.  The 
instrument  used  (ASASP-X)  had  been  subjected  to  laboratory  tests,  first  with 
particles  of  known  sizes  to  determine  accuracy  in  sizing  and  then  with  the 
smoke  of  phosphorus  to  determine  response  to  high  densities.  Calibrations 
with  respect  to  absolute  density  for  each  increment  of  particle  radius  are 
difficult  to  perform  with  accuracy  and  have  not  been  accomplished  for  this 
instrument.  As  stated  then,  the  form  of  the  particle  size  distribution  was 
the  object  of  this  measurement.  TFe 'absolute  magnitude  of  the  absorption  was 
obtained  by  using  another  technique--the  spectrophone. 


optica i 
ven 

...  .  -  ,  o  A 
I  C  *u  wc  u 


Before  the  measurements  Involving  the  spectrophone  are  discussed,  an 
nephelometer  measurement  of  aggregate  particle  scattering  from  within  a  giv 
volume  into  a  fixed  solid  angle  will  be  introduced  since  it  might  be  ,*e!at 
to  the  net  absorption.  The  nephelometer  averages  over  a  volume  as  does  the 
spectrophone;  but,  like  the  particle  measurement  systems  (PMS)  Instrument,  it 
measures  scattering.  The  calibration  and  use  of  a  meteorological  research 
instrument  (MRI)  type  2050B  unit  will  be  discussed  later  in  this  report. 


The  authors'  application  of  spectrophones  to  atmospheric  gases  and  particu¬ 
lates  is  documented  in  the  literature.2  In  a  more  recent  publication,3  the 
authors  discuss  their  spectrophone  research  on  the  smoke  of  white  phosphorus 
(WP).  Spectrophone  theory  has  been  documented,4  5  and  a  variety  of  ASL 


2C .  W.  Bruce  and  R.  G.  Pinnlck,  1977,  "In-situ  measurements  of  aerosol 
absorption  with  a  resonant  CW  laser  spectrophone,"  Ajopl  Opt,  16:1762 

3C.  W.  Bruce  and  Y.  P.  Yee,  1980,  "In-situ  measurement  of  the  ratio  of  aerosol 
absorption  to  extinction  coefficient,"  Appl  Opt,  19:1893 

4E„  L.  Kerr  and  J.  G.  Attwood,  1968,  "The  las  ■  Illuminated  absorptivity 
spectrophone:  A  method  for  measurement  of  weak  absorptivity  in  gases  at  laser 
wavelengths,"  Appl  Opt,  7:915 

5l.  B.  Kreutzer,  1971,  "Ultralow  gas  concentration  infrared  absorption 
spectroscopy,"  Appl  Phys,  42:2934 
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applications  Is  described  In  an  ERADCOM  report.6  The  particular  system  used 
for  these  measurements  was  designed  for  field  use  and  is  tunable  and 
stabilized  for  several  10wm  laser  lines.  Figure  1  shows  a  schematic  cross 
section  of  this  system,  and  figures  2  through  6  show  on-site  photographs  of 
the  instrumentation.  The  Instrumentation  was  first  mounted  within  the  trench, 
but  for  the  later  data,  the  Instrumentation  was  located  outside  and  adjacent 
to  the  polyethylene  cover,  sampling  within  the  trench  through  tubes. 


AMBIENT  ADSORPTION 

The  ambient  gas  and  particulate  absorption  coefficients  were  measured  before 
the  tests.  The  particulate  absorption  was  obtained  Indirectly  by  using  parti¬ 
cle  counting  results.  Lorenz-MIe  calculations  were  based  on  these  data  and 
complex  refractive  indices  for  soil -based  (clay)  dust.  Though  not  really 
precise,  this  technique  usually  Is  accurate  within  a  factor  of  three.  The 
ambient  level  of  particulate  absorption  (coefficient)  was  found  to  be  approxi¬ 
mately  10-3  km-1.  The  total  (gaseous  and  particulate)  absorption  coefficient 
as  measured  by  the  snectrophone  was,  as  expected,  much  higher.  Conditions 
within  the  trench  were  relatively  repeatable  from  day  to  day.  Typical  param¬ 
eters  for  calculations  and  calculated  values  are  shown  In  table  2. 


TABLE  2.  CALCULATED  RANGE  OF  GASEOUS  ABSORPTION  COEFFICIENTS, 
a,  IN  TRENCH  (PRE-SMOKE  TEST  VALUES) 


Date/ 

Dew-Point 

Temperature 

Partial 

Pressure 

H,0  V 

co2 

@330  ppm 

a 

Total 

(minus  trace  gases) 

Time 

(°C) 

(to rr) 

(km-1) 

(km"1) 

(km"1) 

28-29  Jul 

4-5  pm  80 

26.2 

0.66 

0.08 

0.74 

30  Jul 

5  pm 

89.0 

35.0 

1.11 

0.08 

1.19 

3  Aug 

11  am 

88.8 

34.8 

1.09 

0.08 

1.17 

6Y.  P.  Yee,  C.  W.  Bruce,  and  R.  J.  Brewer,  1980,  Gaseous/Parti cul ate  Absorp- 
tlon  Studies  at  WSMR  Using  Laser  Sourced  Spectrophones,  ASl-TR-0065,  US  Army 
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Figure  4.  This  view  (just  above  that  of  figures  2  and  3)  shows  the  MR  I  nephelomet 
and  PM 3  counter  and  dew-point  hygrometer  electronics  in  the  foreground. 


Spectrophone  calibration  setup  in  laboratory  near  site. 


The  partial  wa ter  pressure  was  based  on  extrapolations  front  the  data  of  Yin  et 
al,7  and  the  results  were  based  on  data  of  Shumate  et  al ,B  about  11  percent 
lower  for  these  partial  pressures.  0;te  tunnel  section  (within  4  m  of  the 
measurements)  was  removed  for  the  28-29  August,  4  to  5  pm  times.  The  section 
cover  was  reinstalled  before  the  aerosol  tests,  and  the  data  here  for  30  July 
and  3  August  are  for  completely  covered  trenches. 

The  measured  gaseous  absorption  coefficient  corresponding  to  the  first  example 
(data  of  28-29  July  1980)  was  1.1  km-1.  Probable  error  figures  on  this  mea¬ 
surement  are  not  precise  but  are  estimated  to  be  0.2,  largely  because  of 
temporal  variations  due  to  unknown  source(s).  Higher  than  predicted  values  of 
the  absorption  coefficient  frequently  occur  In  ASL  field  measurements  of 
ambient  absorption  at  lOum  wavelengths. 

Absorption  due  to  trace  constituents  could  elevate  the  value  considerably.  A 
rough  rule  of  thumb  Is  that  trace  gas  absorption  (within  an  absorption  band 
region)  often  reaches  or  exceeds  1  (ppm-km)-1  and,  of  course,  1  ppm  Is  not  a 
high  concentration. 

During  the  measurement  program,  test  project  personnel  frequently  mentioned 
the  desire  to  obtain  total  absorption  values  lower  than  those  measured  ojr_ 
predicted  (minus  trace  gas  contributions).  Another  reason  for  Interest  In  the 
ambient  gaseous  absorption  Is  that  the  smoke  aerosol  level  approached  these 
levels  near  the  end  of  each  test. 


AEROSOL  CHARACTERIZATION  FOR  TESTS 

Aerosol  characterization  for  two  test  situations  will  be  discussed.  The  dates 
and  approximate  starting  times  are  listed  as  the  latter  two  examples  in  table 
2  (Information  for  30  July  1980  and  3  August  1980).  The  most  notable  dif¬ 
ference  In  the  conduct  of  these  two  tests  Is  the  time  of  day.  Aerosol  char¬ 
acteristics  were  similar.  In  th-  mean. 

The  procedure  was  to  Ignite  phosphorous  smoke  grenades  In  the  trench  at  about 
the  one-quarter  and  three-quarter  length  positions.  The  smoke  dispersed 
roughly  over  a  2-h  period  during  which  aerosol  measurements  were  made.  The 
data  and  analyses  will  be  presented  as  functions  of  the  evolving  aerosol 
system:  and  Implications  for  nonlinear  effects  will  be  discussed. 

Particle  size  distribution  evolved  generally  as  expected  from  prior  laboratory 
measurements  on  the  same  aerosol;  that  Is,  the  mean  densities  decreased  and 
the  peak  moved  to  progressively  smaller  sizes.  Figure  7  shows  a  composite 
graph  of  the  distribution  pattern  for  the  first  test.  These  data  represent 
roughly  1-mln  averages  (data  were  not  always  continuously  available  In 
time).  Circulation  under  the  loosely  fit  trench  cover  (and  pumping  action  by 


7P.  K.  L.  Yin  and  R.  K.  Long,  1968,  “Atmospheric  absorption  at  the  line  center 
of  P ( 20 )  C02  laser  radiation,"  Appl  Opt  7:1551 

®M.  S.  Shumate  et  al ,  1976,  "Water  vapor  absorption  of  carbon  dioxide  laser 
radiation,"  Appl  Opt  15:2480 
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fluctuations  In  air  pressure  on  top  of  the  plastic  cover)  was  sufficient  that 
any  type  of  settling  calculations  to  determine  decreasing  densities  would  be 
questionable. 

The  differential  absorption  coefficients  (absorption  per  micrometer  particle 
radius)  were  calculated  by  using  '.orenz-MIe  theory  with  optical  constants 
appropriate  to  the  composition  of  the  smoke.  To  accomplish  this  calculation, 
the  relative  humidity  for  the  test  was  computed  from  the  measurements  of 
figure  8.  The  relative  humidity  was  used  In  turn  to  calculate  the  percent  of 
water  for  the  hygroscopic  smoke  on  the  basis  of  measurements  made  In  our 
laboratory  and  also  predicted  by  the  theory  of  Hanel  and  Bull  rich.9  The 
relationship  Is  plotted  In  figure  9.  finally,  the  complex  indices  were 
obtained  from  a  work  by  Querry  and  Tyler,10  but  heir  data  posed  a  problem; 
two  sets  of  complex  indices  were  obtained  from  _.<1s  source  (figure  10).  We 
have  discussed  the  strange  appearing  results  with  the  authors  who  maintain 
that  both  sets  are  correct.  We  have  simply  used  them  (no  measurements  fell  In 
the  disjointed  data  gap). 

Figure  11  represents  a  plot  of  the  differential  absorption  based  on  the  data 
of  figure  7.  The  absorption  peak  occurs  at  a  radius  of  about  1.2um  soon  after 
the  smoke  Is  dispersed  and  decreases  relatively  smoothly  to  about  0.4ym  in 
roughly  2  h;  in  the  same  time  the  peak  absorption  level  has  dropped  roughly 
two  orders  of  magnitude.  At  the  midpoint  In  this  2-h  evolution,  80  percent  of 
the  absorption  occurs  within  a  radius  span  of  approximately  0.30um. 

The  background  particulate  absorption  Is  represented  by  the  lowest  curve 
labeled  1641,  that  Is,  before  the  smoke  was  released.  The  contribution  here 
(at  about  10"3km-1)  Is  well  below  the  ambient  gaseous  absorption. 

The  total  absorption  coefficient  as  calculated  from  the  particle  size  and 
density  Information  is  plotted  as  a  function  of  time  In  figure  12.  For  this 
presentation,  i-min  averages  obtained  at  approximately  10-mln  intervals  are 
connected  by  straight-line  segments.  These  values  are  much  lower  than  those 
of  the  spectrophone  measurements.  The  spectrophone  yielded  absorption  and 
extinction  coefficient  Deak  values  (In  separate  measurements)  of  about  200  and 
220  km"1,  while  the  par  dele  counter  result  was  about  10  km"1.  This  differ¬ 
ence  Is  believed  to  be  due  primarily  to  low  counting  efficiencies  and  sampling 
errors  In  unknown  proportions.  As  the  purpose  of  sizing  and  counting  for 
these  tests  was  to  obtain  measurements  of  the  form  of  the  size  distribution  as 
a  function  of  time  rather  than  absolute  values,  this  difference  Is  not  a  prob¬ 
lem.  Prior  calibration  of  the  counter  using  monodisperse  aerosols  precludes 
sizing  as  a  prime  source  of  error.  Likewise,  the  optical  constants  and  calcu- 
latlonal  scheme  introduce  uncertainties  that  are  very  small  compared  with  the 
difference. 


9G.  Hanel  and  K.  Bullrlch,  1978,  "Physico-chemical  property  models  of 
tropospheric  aerosol  particles,"  Beltrage  Zur  Physlk  der  Atmosphare,  51:129 

1UM.  R.  Querry  and  I.  L.  Tyler,  1978,  "Complex  refractive  Indices  in  the 
infrared  for  H.JP04  In  water,"  J  Opt  Soc  Am,  68:1404 


18 


140 


time  of  day 

Figure  8.  Some  temperatures  and  dew-point  temperatures  obtained 
during  the  first  test. 


Figure  9.  Growth  of  particles  with  relative  humidity  (smoke  of 
phosphorus)  O.lym  to  10pm  (radius)  particles. 


Figure  10.  Plot  of  measured  complex  indices  for  phosphoric  acid 
at  10.6pmx. 
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differential  absorption  coefficient  ik« 


Figure  11.  Differential  absorption  for  the  first  test 


A  comparison  of  figures  7  and  11  shows  that  the  largest  size  particles  are 
most  important  in  the  calculation  of  the  total  absorption  because  of  the  rapid 
increase  with  particle  radius  of  the  Lorenz-Mie  efficiency  factor  for  the 
absorption  cross  section.  Therefore,  the  apparent  Increase  with  time  in  the 
densities  of  the  relatively  small  particles  (whose  signals  are  probably  masked 
in  the  counter  by  those  from  the  larger  particles)  is  not  of  concern  here.  A 
designed-in  characteristic  of  the  counter  is  also  partly  at  fault  in  this 
misrepresentation  which  causes  a  "tailing-off"  effect  in  the  signals  repre¬ 
senting  progressively  smaller  sizes. 

Calculations  to  follow  are  based  on  the  peak  absorption  values  which  aqree 
well  with  the  author's  predictions  based  on  the  laboratory  measurements  of 
reference  3  (within  15  percent). 

The  measurement  of  extinction  coefficient  derived  from  power  less  in  the 
spectrophone  is  satisfactory  at  high  values  (early  in  the  test)  but  is  not 
useful  below  10  km-1  due  to  the  question  of  power  meter  drift.  Figure  13 
snows  the  extinction  coefficient  for  the  test  of  figure  12.  The  slope  of  the 
decreasing  extinction  coefficient  with  time,  however,  agrees  well  with  that 
found  by  using  the  PMS  particle  counter  data.  An  electronic  problem,  presumed 
to  be  overheating  of  system  electronics  prevents  complete  comparison  of  spec¬ 
trophone  absorption  data  for  this  test. 

Time  variations  culled  from  the  extinction,  the  density  distribution,  and  the 
absorption  measurements  at  the  location  of  the  apparatus  were  greater  than  a 
factor  of  five  for  Integration  times  less  than  3  s  and  about  a  factor  2  to  3 
for  Integration  times  of  10  s.  However,  a  spatial  density  average  over  as 
much  as  10  percent  or  greater  volume  of  the  trench  probably  varies  less  than 
20  percent. 

The  temperature  time  series  for  the  second  test  is  given  by  figure  14. 

Temperatures  in  the  trench  were  somewhat  lower  than  the  temperatures  in  the 
first  test  even  though  the  second  tests  were  during  midday  (the  30  July  test 
was  later  in  the  afternoon).  Particle  size  distributions  and  their  evolution 
in  time  were  similar  in  form  though  the  corrosive  action  of  the  smoke  caused  a 
reduction  in  particle  counter  laser  power  for  this  test.  The  spectrophone 
measured  peak  absorption  and  extinction  values  were  virtually  the  same  at 
200  km-1.  Absorption  coefficient  as  a  function  of  time  for  this  test  is  shown 
in  figure  15.  For  a  resolution  time  of  about  2  s,  the  magnitude  fluctuates  by 
more  than  an  order  of  magnitude  (figure  16).  The  same  data,  when  averaged 
over  about  a  minute  for  each  of  the  sampling  periods,  yield  a  much  more  steady 
decrease.  This  decrease  has  nearly  the  same  mean  slope  as  that  calculated 
from  the  particle  counter  data  though  the  magnitudes  are  very  different. 

In  contrast  to  the  spectrophone,  the  nephel ometer  measurement  yields  primarily 
particle  scattering  information  (molecular  scattering  is  generally  much 
lower).  Due  to  their  high  sensitivity,  data  were  obtained  only  at  the  very 
end  of  these  tests--when  the  ends  of  the  trench  were  opened  near  test  comple¬ 
tion.  The  nephelometer  output  (scattering  c  'fficient  at  visible  wavelengths) 
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TIME  OF  DAY 


Figure  15.  Absorption  coefficient  as  a  function  of  time  for  the 
second  test. 


Figure  16.  Illustration  of  large  fluctuations  in  the  absorption 
coefficient  (2  s  time  resolution)  early  during  the 
second  test. 


?A 


was  related  in  a  laboratory  calibration  to  extinction  and  absorption  coeffi¬ 
cients  at  a  wavelength  of  10.6wn  (figure  17).  Since  there  was  no  evidence  of 
nonlinearity,  a  single  constant  serves  to  relate  the  volume  scattering  coeffi¬ 


cient  to  extinction,  that  is,  i-  — — -  and  another  to  relate  the  scat- 

Sc.  Coef. 


tering  coefficient  to  absorption  0.96  x  -i-  — — -  for  the  smoke  of  WP. 

L  F  Sc.  Coef . J 

Figure  18  shows  how  the  particulate  concentration  can  fall  at  the  center 
position  when  the  wind  is  allowed  to  flush  the  trench.  A  nephelometer  near 
the  target  end  of  the  trench  (about  6  m  from  the  end),  though,  recorded  that 
the  fluctuating  smoke  attenuation  decreased  in  that  region  rapidly  to  below 
2  km'1  at  about  1218,  while  near  the  center,  the  absorption  coefficient  level 
was  still  near  10  km"1.  For  about  1  min  during  which  laser  pulse  7718  was 
fired  through  the  trench  (1214  to  12i„),  this  nephelometer  measured 
a  =  1.51  ±  0.C6  km"1. 

The  smoke  of  WP  is  hygroscopic,  quickly  drawing  water  from  the  vapor  state. 
The  reduction  of  the  relative  humidity  resulting  from  the  introduction  of  the 
smoke  into  the  tunnel  was  calculated  on  the  basis  of  prior  laboratory  measure¬ 
ments,  that  is,  a  smoke  production  efficiency  actor  of  kE  =  0.35  m2/gm  and  a 
fractional  growth  in  particle  radius  of  r/rQ  1.21  for  an  initial  relative 
humidity  of  35  percent.9  By  using  the  initial  value  of  absorption  to  calcu¬ 
late  the  mass  density  of  the  smoke  (pm  =  0.16  gm/m3)  and  assuming  a  dew-point 
temperature  of  84°F  to  obtain  the  initial  partial  pressure  of  water  (29.8 
torr),  a  partial  pressure  change  of  less  than  0.1  torr  was  obtained.  This 
change  is  negligible  for  effects  of  interest  here. 


EFFECT  OF  EVAPORATIVE  CLEARING  ON  HIGH  POWER  BEAMS 

The  tendency  to  punch-through  the  absorbing  countermeasure  aerosol  WP  is  of 
prime  interest  here,  and  parameters  from  the  two  similar  tests  will  be  used  as 
a  basis.  The  thermodynamic  and  optical  bases  for  the  calculations  will  be 
applied  by  using  laser  beam  parameters  relevant  to  a  system  whose  beam  param¬ 
eters  are  to  be  described  (US  Army  Missile  Command  [MICOM]  S3  system). 


9G.  Hanel  and  K.  Bull  rich,  1978,  "Physico-chemical  property  models  of 
tropospheric  aerosol  particles,"  BeUrage  Zur  Physik  der  Atmosphare,  51:129 
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Figure  17.  Settling  chamber  measurements:  nephelometer  output 
relationship  with  extinction/absorption  coefficients. 
Absorption  coefficient  =  0.96  x  extinction  coefficient. 


TIME  OF  DAY 

Figure  18.  Absorption  coefficient  at  the  end  of  the  second  test, 

first  by  the  center  spectrophone  and  then  by  the  center 
nephelometer. 
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The  thermodynamical  calculations  for  the  clearing  effect  begin  with  the  dic¬ 
tionary  of  terms  listed  In  table  3.  Glickler11  and  Sutton*2  used  elaborate 
(but  still  Inexact)  integral  solutions  to  describe  clearing  by  fo£. 


TABLE  3.  DEFINITION  OF  SYMBOLS  FOR  EVAPORATIVE  CLEARING 


Q 

rp 

rb 

r 

P 

E 


Pm 

c 


H 


i 

t. 


heat  (joules) 

radius  of  particle  (micrometers,  centimeters) 
radius  of  beam  (centimeters) 
variable  representing  radius 
instantaneous  power  (watts) 
beam  energy  (joules) 

bulk  density  of  particulate  material  (grams/cubic  centimeters) 

specific  heat  of  particulate  material  ( joules/grams  degrees 
Kelvin) 

latent  heat  of  vaporization  (joules/grams) 

Lorenz-Mie  efficiency  factors  for  absorption  and 
extinction 

mean  temperature  of  particle  (degrees  Kelvin) 
time  (seconds) 


Note  that  pm,  c,  Tg0^,  H,  and  the  complex  refractive  indices  leading  to  Qa 
and  Qe  are  functions  of  the  relative  humidity.  The  values  are  listed  in 
appendix  A. 


ll$.  L.  dicker,  1971,  "Propagation  of  a  10. 6pm  laser  through  a  cloud  includ¬ 
ing  droplet  vaporization,"  Appl  Opt,  10:644 

12G.  W.  Sutton,  1978,  "Fog  hole  boring  with  pulsed  high  energy  lasers:  An 
exact  solution  including  scattering  and  absorption,"  Appl  Opt,  17:3424 


27 


CALCULATIONS,  PART  I:  EVAPORATION  OF  SINGLE  PARTICLES 


Particles  are  first  heated  to  boiling  T,  assuming  spatially  uniform  heating  of 
the  particle,  and  then  energy  Is  added  to  evaporate  them. 

To  reach  boiling  temperature,  the  rate  of  heating. 


r  p 

Q  =  nR  2Q - beam  power  density 

p  °L("rb2) 


/R  \2  t 

then  is  Integrated,  Q  =yjpj  Qa  j  Pdt,  and  equated  with  the  energy 


required  to  boll  the  particle. 


Q  =  T  nRp3  Dm  rAT  ‘ 

Now  j  Pdt  is  just  the  energy,  E,  which  has  passed  the  plane  of  the 
o 

particle  considered  (to  time  t  ).  Solving  for  this  energy,  one  obtains 


Eg  .  (**)  AT  . 

Q 


To  add  the  heat  of  vaporization, 


c  _  r4n->  pm  RpRB2  „ 

tLHV  "  - - - H  * 

a 


The  sum  of  these  two  processes  represents  the  total  beam  energy  which  has 
passed  the  plane  of  a  particle  of  radius  Rn  at  the  time  of  evaporation,  that 
is,  p 


WV  ■  MV  *  EEHv(Rp) 


An  examination  of  the  S3  laser  pulse  intensity  profiles  led  to  the  suggestion 
that  the  pulse  be  modeled  as  triangular  in  time.*  Samples  pulse  shapes  are 
shown  in  appendix  B.  Figure  19  shows  sample  real  and  approximate  pulse  forms. 


*Bill  Jones,  MICOM,  private  communication 
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Figure  19.  Sample  real  and  approximate  pulse  forms. 
Sample  (and  typical)  parameters  are 


E( total )  =  i°ules 

‘pulse  =  3'2l,s 
fmax  ‘  125 

The  objective  here  is  to  determine  the  time  to  evaporate  particles  and  then  to 
calculate  the  energy  remaining  for  each  power  lobe  of  the  far  field  pattern. 

t' 

E  =  f  Pdt  . 
cum  J 


Then  using  a  linear  form 


*CttW  *  /  ("a  + 


=  2-  (t’)2  +  Pra,vt'  where  m  - 

l  ITlaX 


pmax 

Sulse 
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Solving  this  for  t* ,  the  cumulative  pulse  time. 


Here  It  Is  assumed  that  the  particle  Is  vaporized  and  vaporous  at  time  t* . 
The  first  simplifying  assumption  1?  justified  since  clearing  Is  effective  when 
particle  diameters  are  reduced  In  boiling  by  considerably  less  than  an  order 
of  magnitude  (see  figure  11).  The  second  assumption  follows  from  kinetic 
calculations  showing  that  dispersal  velocities  are  sufficient  to  form  an 
effective  vapor  state  within  an  Interval  very  short  compared  with  the  pulse 
duration. 

Figure  20  shows  the  problem  schematically. 


t 


TIME  TO  REACH 
FOILING  TEMP. 


-  ENERGY  REMAINING 
TO  PUNCH -THRU 
(OUR  INTEREST  HERE) 


ADDING  LHV 


Figure  20.  Expenditure  of  energy. 


Now  the  energy  (or  power)  density  distribution  in  the  beam  must  be  determined. 

The  job  is  simplified  since  the  laser  output  is  a  line-square,  that  Is,  two 
pairs  of  parallel  slits  in  orthogonal  direc cions  (figure  21).  A  sample  burn 
pattern  (quite  overburned)  Is  shown  in  appeidix  C. 


Figure  21.  Near,  near  field  pattern. 


The  power  densities  In  the  near,  near  field  are  therefore  those  of  orthogonal 
pairs  of  finite  slits  and,  In  the  far  field,  are  the  combined  diffraction 
patterns  of  these  slits.  Actual  burn  patterns  (sample  In  appendix  D)  look 
quite  like  the  sketch  of  figure  22. 
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(There  will  be  some  energy 
C,  in  each  quadrant,  but  the 
amount  will  be  relatively 
small  compared  with  the 
identified  energy  "lobes.") 


Figure  22.  Far  field  pattern. 


The  approach  Is  to  estimate  the  distribution  of  energy  In  this  pattern  or  that 
In  the  area  under  the  twin  slit  Intensity  profile. 


I 


t 


that  is,  the  single  slit  pattern  (first  term)  modulated  by  that  of  the  twin 
slits.  Beta  and  alpha  contain  the  viewing  angle  and  dimensional  parameters  of 
the  calculation. 


Figure  23  shows  a  perspective  sketch  of  the  far  field  beam  profile. 
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P  -  KEl^a^  -  KaE^I^  , 


where  Ij  are  the  relative  Intensities  and  the  K  is  found  from  the  total  energy 
and  the  pulse  shape. 

Appendix  E  contains  the  details  of  the  calculation  of  the  distribution  of 
energy  In  the  laser  beam  to  which  this  formula  Is  applied. 

The  distribution  of  energy  In  the  lobes  as  calculated  in  appendix  E  Is  given 
in  table  4. 


TABLE  4.  DISTRIBUTION  OF  ENERGY  IN  FAR  FIELD  BEAM  LOBES  FOR  S3  LASER 


Lobe* 

Number 

Percent  E 
per  Lobe 

Percent  E  per  Lobe 
Lobe  Area 
(cm-2) 

Energy  Density  per  Lobe 

for  2 total  ■  200  J 
(J/cm2) 

0 

18.7 

7.19 

14.4 

1 

8.5 

4.47 

8.94 

2 

6.3 

3.32 

6.64 

3 

3.7 

1.95 

3.90 

4 

1.5 

0.79 

1.58 

5 

0.3 

0.16 

0.32 

♦Increasing  away  from  center 


The  total  area  Involved  at  the  target  distance  Is  approximately  40.6  cm2. 

The  clearing  time  and  then  the  fractional  energy  remaining  to  be  transmitted 
after  clearing  will  be  calculated.  The  clearing  time  and  fractional  energy 
will  then  be  obtained  for  the  entire  beam. 

The  clearing  process  effectively  reduced  the  absorption  to  gaseous  levels 
which  are  much  lower,  that  Is,  to  near  normal  atmospheric  levels.  The  scat¬ 
tering  is,  of  course,  also  greatly  reduced;  however,  for  the  smoke  of  WP,  the 
fraction  of  C02  laser  beam  energy  scattered  to  that  absorbed  Is  less  than  10 
percent.  The  atmospheric  propagation  Is  therefore  related  to  the  punch- 
through  effect  In  a  form  Illustrated  In  figure  24. 
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Figure  24.  Illustration  of  relationship  between  laser 
pulse  and  propagation  parameters. 


Tables  5  and  6  show  rapid  and  effective  clearing  as  a  function  of  smoke  parti¬ 
cle  size  for  a  relative  humidity  of  35  percent.  Earlier  data  show  that  very 
few  particles  are  larger  than  Rp  :  4um  and  that  the  clearing  efficiency  is' 

almost  independent  of  size  to  that  point  because  of  the  increase  in  Lorenz-Mie 

efficiency  with  radius  (up  to  Rp  a  4um)  in  spite  of  the  growing  mass  to  be 

evaporated.  Coincidentally,  this  convenient  independence  disappears  for  the 
larger  particles.  The  last  column  represents  the  efficiency  of  punch-through 
or,  in  other  words  the  percentage  of  energy  in  the  cleared  field.  The  trans¬ 
mission  Is  not  unity  at  that  point  {as  was  mentioned)  but  reverts  nearly  to 
the  pretest  levels. 
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TABLE  5.  PARAMETERS  OF  SINGLE  PARTICLE  CLEARING  AS  FUNCTIONS  OF  PARTICLE 
SIZE  AND  BEAM  LOBE  NUMBER  FOR  ETOm  =  200  J  AND  Tpu,  «  =  3.2Ms. 


Center  Lobe  (0) 
Energy  =  40.0  J 


<k> 

Qa 

Time  to 
Evaporate 
(ys) 

0.01 

0.00202 

1.68  x  10"1 

0.05 

0.0101 

1.68  x  10“* 

0.10 

0.0203 

1.67  x  10"1 

0.15 

0.0305 

1.67  x  10"1 

0.20 

0.0511 

1.32  x  10“i 

0.40 

0.0828 

1.64  x  10_i 

0.60 

0.127 

1.60  x  lO'i 

1.0 

0.222 

1.52  x  10"i 

2.0 

0.490 

1.38  x  10‘i 

4.0 

0.904 

1.50  x  10"! 

8.0 

1.217 

2.25  x  10"i 

10.0 

1.265 

2.72  x  lO'i 

13.0 

1 . 2600 

4.20  x  10“i 

20.0 

1.2123 

6.00  x  10"i 

25.0 

1.1663 

8.09  x  lO'i 

30  0 

1.1663 

1.01 

40.0 

1.1663 

1.47 

50.0 

1.1663 

2.11 

60.0 

1.1663 

>3.2 

Fraction  of  Energy  j 

In  Remaining  i 

Portion  of  Pulse  j 

0.90  j 

j| 

0.90  |j 

0.90  ; 

0.90  t 

0.92  j 

0.90  \ 

0.90  | 

0.91  ; 

0.92 

0.91 

0.87 

0.84 

\ 

0.76 

0.66  | 

0.56  ) 

l 

0.47 

0.29 

0.12 

0.00 
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TABLE  5.  (Cont) 


First  Side  Lobe  (l) 

Energy  =  18.1  J 

A) 

Time  to 
Evaporate 
(ys) 

Fraction  of  Energy 
In  Remaining 
Portion  of  Pulse 

0.01 

2.76 

X 

10-1 

0.83 

0.05 

2.76 

X 

10-1 

0.83 

0.10 

2.75 

X 

10-1 

0.84 

0.15 

2.74 

X 

10“  1 

0.84 

0.20 

2.16 

t 

X 

10-1 

0.87 

0.40 

2.69 

X 

10“1 

0.84 

0.60 

2.63 

X 

10"! 

0.84 

1.00 

2.50 

X 

10"  i 

0.0 

2.00 

2.26 

X 

lO'i 

4.00 

2.46 

X 

io-i 

0.0 

Second  Side  Lobe  (?) 

Energy  =  13.5  J 

0.01 

3.77 

V 

lO'i 

0.78 

0.05 

3.77 

X 

lO'i 

0.78 

0.10 

3.77 

X 

10"i 

0.78 

0.15 

3.77 

X 

lO'i 

0.78 

0.20 

2.94 

X 

io-i 

0.82 

0.40 

3.67 

X 

lO'i 

0.0 

0.60 

3.58 

X 

io-i 

0.0 

1.00 

3.41 

X 

10"! 

0.0 

2.00 

3.07 

X 

lO'i 

0.0 

4.00 

3.34 

X 

10*1 

0.0 

TABLE  5.  (Cont) 


Third  Side  Lobe  (1) 
Energy  =  7.88  J 


1% 

Time  to 
Evaporate 
( us) 

Fraction  of  Energy 
In  Remaining 
Portion  of  Pulse 

0.01 

6.75  x  lO’1 

0.62 

0.05 

6.75  x  10'1 

0.62 

0.10 

6.75  x  10'1 

0.62 

0.15 

6.75  x  10'1 

0.62 

0.20 

5.20  x  10”*1 

0.70 

0.40 

6.57  x  10'1 

0.0 

0.60 

6.41  x  10"1 

0.0 

1.00 

6.07  x  10'1 

0.0 

2.0  0 

5.44  x  10'1 

0.0 

4.00 

5.95  x  nr1 

0.0 

Fourth  Side  Lobe 

Energy  =  3.20  J 

0.01 

2.37 

0.69  x  10"i 

0.05 

2.37 

0.69  x  10'1 

0.10 

2.37 

0.69  x  10'1 

0.15 

2.37 

0.69  x  lO"1 

0.20 

1.56 

0.26 

0.40 

2.04 

0.91  x  10'1 

0.60 

2.14 

0.0 

1.00 

1.95 

0.0 

2.00 

1.66 

0.0 

4.00 

1.89 

0.0 

(Particles  In  the 

fifth  side  lobe  require  >  3.2vs) 
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: 

AND  USING  PARTICLE  RADII  OF  PEAK 

ABSORPTION 

k-  i 

i 

Time:  1159 

Run  Number  7716 

< 

i 

“  \ 

[• 

w  ■  3-0l,s 

Etotal  =  144-9  J 

Rp  =  l.Oum 

Qa  =  0.2217 

, 

k .  4 

Beam  Lobe 

Time  to  Boil 

Fraction  of  Energy 

Number 

Pmax  ancl  EvaPorate 

(W)  (ms) 

in  Remaining  j 

*  1 

r* 

(from  center) 

Portion  of  Pulse  j 

4 

- i 

TABLE  6.  (Cont) 

Time:  1214 
Run  Number  7718 
tmax  =  3.3ps 

Etotal  =  204,5  J 
Rp  =  0.7pm 

Qa  =  0.1469 


Beam  Lobe 
Number 

(from  center) 

>wV 

Time  to  Boil 
and  Evaporate 
(ps) 

Fraction  of  Energy 
in  Remaining 
Portion  of  Pulse 

0 

2.32  x 

10' 

1.73  x  10'1 

0.90 

1 

1.05  x 

107 

2.85  x  10-1 

0.84 

2 

7.81  x 

106 

3.89  x  10"i 

0.78 

3 

4.59  x 

106 

6.96  x  10"i 

0.62 

4 

1.86  x 

106 

2.44 

0.65  x  10"i 

5 

3.72  x 

105 

>3.3 

0.0 

Time:  1221 
Run  Number  7719 

^ax  =  3,0)jS 

Etotal  =  200  J 
Rp  =  0.6pm 

Q  =  0.1274 


0 

2.49 

x 

107 

1.59 

X 

10"i 

0.90 

1 

1.13 

X 

107 

2.61 

X 

lO'i 

0.83 

2 

8.40 

X 

106 

3.58 

X 

lO'i 

0.78 

3 

4.93 

X 

106 

C\J 

• 

vo 

X 

io-i 

0.62 

4 

2.00 

X 

10b 

2.28 

0.66  x  10"  1 

5 

4.00 

X 

105 

>3.0 

0.0 
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TABLE  6.  (Cont) 


Time:  1227 
Run  Number  7720 

Wx  ”  3.3ns 

Etota1  =  198  J 
Rp  =  0.45um 

Qa  =  0.0937 


Beam  Lobe 
Number 

(from  center) 

PiW 

Time  to  Boil 
and  Evaporate 
(us) 

Fraction  of  Energy 
in  Remaining 
Portion  of  Pulse 

0 

2.24  x 

107 

1.81  x  10" 1 

0.89 

1 

1.02  x 

107 

2.96  x  10"1 

0.83 

2 

7.56  x 

106 

4.06  x  10"1 

0.77 

3 

4.44  x 

10s 

7.30  x  10"1 

0.60 

4 

1.80  x 

106 

2.73 

0.34  x  10"1 

5 

3.60  x 

105 

>3.3 

0.0 

Time:  1232 
Run  Number  7721 

^max  =  3,7wS 

Etotal  “  193 *2  J 
Rp  =  0.40ym 

Qa  =  0.0828 


0 

1.95 

X 

107 

2.09 

X 

10"1 

0.89 

1 

8.88 

X 

106 

3.42 

X 

10"1 

0.82 

2 

6.58 

X 

IQ6 

4.70 

X 

10" 1 

0.76 

3 

3.86 

X 

106 

8.48 

X 

10" 1 

0.59 

4 

1.37 

X 

10'6 

3.53 

0.34  x  10*2 

5 

3.13 

X 

105 

>4.0 

0.0 
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CALCULATIONS,  PART  II:  COLLECTIVE  EFFECT 


From  these  results  one  can  calculate  the  net  clearing  effect  on  the  beam. 
This  punch-through  phenomenon  will  be  illustrated  in  the  following  calcula¬ 
tions.  For  this  purpose  the  particle  size  corresponding  to  the  peak  absorp¬ 
tion  for  selected  settling  times  will  be  used. 

The  expenditure  of  energy  as  a  function  of  aerosol  loading  can  be  calculated 
from 


E(L)  =  E(t  =  0  to  t')e"ax  +  Eft  =  to  tp(jlse) 


P(t)dt  + 


P(t)dt  , 


where  t'  is  the  time  to  evaporate  and  t_uise  is  the  total  pulse  duration. 
Here  t'  will  have  been  increased  over  that  of  the  single  particle  function  as 
a  result  of  the  mass  loading.  For  particles  of  given  size  and  composition,  t' 
is  an  inverse  function  of  Pmax-  (L)  represents  the  path  length  within  the 
aerosol  medium.  The  time  to  evaporate  at  a  given  position  in  the  measurement 
is  now  dependent  on  that  position  although  the  above  equation  is  written  as 
though  ic  depends  only  on  the  end  point.  This  simp! ificatloii  does  not  signif¬ 
icantly  alter  the  results  here.  The  smoke  of  phosphorus  clears  significantly 
for  a  fairly  broad  range  of  aerosol  absorption  coefficients  (or  alternatively 
a  range  of  aerosol  mass  loading  values).  For  relatively  low  absorption  coef¬ 
ficients,  clearing  can  do  little  to  enhance  the  total  energy  transmitted.  The 
enhancement  due  to  clearing  will  be  determined  by  the  above  presentation  and, 
using  this  technique,  can  be  extended  to  absorption  coefficients  of  several 
tens  per  kilometer  (typical  field  values)  for  even  the  relatively  low  energy 
system  of  these  tests.  Clearing  continues  to  be  significant  for  this  (and 
other  volatile)  aerosols  to,  roughly,  the  loading  at  which  beam  energy  is 
consumed  in  the  evaporation  process.  Screening  of  available  energy  for  parti¬ 
cles  farther  from  the  source  causes  an  Increase  in  t'  at  about  25  km-1  (effec¬ 
tively)  though  the  increase  depends  somewhat  on  factors  such  as  particle  size 
and  relative  humidity.  For  the  calculations  to  follow,  open  air  dispersal  of 
WP  was  assumed;  therefore,  the  peak  absorption  Is  for  particles  of  inn  to 
1.2 urn  radius  (see  early-time  data  of  figure  11).  For  all  wests  In  the  trench, 
particle  radii  became  progressively  smaller  at  similar  rates;  and  the  attenua¬ 
tion  also  decreased  (on  the  average)  throughout  the  test.  The  incorporatln  of 
the  bulk  attenuation  into  the  calculation  of  t'  yields  the  results  of  figure 
25.  Table  5  shows  that  particle  size  and  Lorenz-Mie  efficiencies  vary  such  as 
to  hold  t',  the  time  to  evaporate  the  particles,  nearly  constant  (independent 
of  particle  size).  Of  course,  particles  with  large  radii  require  more  energy 
to  evaporate.  The  form  of  the  curve,  that  is,  curving  over  quickly  at  the  top 
(on  a  semi log  plot)  augurs  for  the  efficiency  of  the  punch-through  effect 
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Figure  25.  Time  to  evaporate  center  lobe  in  loaded  aerosol  medium. 

Loading  given  by  absorption  coefficient.  Etota-|  =  40.0  J 
(in  lobe),  Rp  =  lum,  tnuise  =•  3.2us. 


since  It  conveys  that  the  time  to  evaporate  increases  rapidly  with  the  absorp¬ 
tion  coefficient  only  for  relatively  high  values  of  the  coefficient.  The 
remaining  step  is  to  show  the  relative  Improvement  in  the  energy  transmitted 
when  punch-through  occurs  (table  7  and  figure  26).  The  Improvement  in  energy 
transmitted  Is  truly  sign1ficant--even  though  the  pulsed  laser  on  which  the 
calculated  parameters  were  based  on  a  relatively  small  one. 

Generally,  thermal  blooming  and  evaporative  clearing  effects  are  interactive; 
that  Is,  the  power  density  Is  affected  by  the  thermal  blooming  which  in  turn 
affects  the  evaporative  clearing.  Thermal  blooming  depends  on  beam  parameters 
(power  density,  pulse  length,  and  pulse  shape)  and  atmospheric  parameters 
(gaseous  and  particulate  absorption,  crosswinds,  and  turbulence).  Appendix  F 
shows  that  blooming  effects  on  the  test  beam  are  not  expected  to  be  strong; 
therefore,  for  this  analysis,  these  effects  have  been  ignored. 
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Figure  26.  Total  energy  transmitted  in  center  spot  as  a  function  of 
absorption  coefficient  for  the  WP-smoke  aerosol.  No 
gaseous  absorption  effects  included. 


CONCLUSIONS 

Aerosol  characterization  performed  on  the  smoke  of  WP  was  used  to  calculate 
punch-through  or  evaporative  clearing  on  a  pu'ised  C02  laser  beam.  This  effect 
can  be  significant  for  HEL  use.  High  energy  pulsed  C02  laser  beams  can  effec¬ 
tively  penetrate  red  phosphorus  (RP),  WP  countermeasure  smoke,  though  the 
degree  of  efficiency  depends  on  beam  and  aerosol  parameters.  This  effective¬ 
ness  is  related  to  the  volatility  and  absorbing  qualities  of  the  smoke  parti¬ 
cles.  Repetitive  and/or  long  pulses  can  encounter  strong  thermal  blooming  due 
to  the  high  level  of  aerosol  absorption  which  counteracts  the  clearing 
effect.  This  report  discusses  the  development  of  an  analytical  model  and  its 
application  to  aerosol  parameters  that  derive  from  measurements  in  an  extended 
enclosure  at  the  MICOM  in  Huntsville,  AL,  and  then  to  aerosol  parameters  that 
are  more  typical  of  atmospheric  disnersion  of  the  smoke  (represented  primarily 
by  relatively  large  diameter  particles).  The  laser  beam  parameters  in  both 
cases  are  relevant  to  a  high  energy  research  laser  owned  by  MICOM  and  used  for 
propagation  studies  in  the  above-mentioned  enclosure  during  the  in  situ  aero¬ 
sol  measurments  presented  in  this  report.  Aerosol  absorption  measurements 
used  in  situ  aerosol  spectrophones  developed  at  ASL.  Thermal  blooming  (the 
stated  original  objective  of  the  MICOM  exploratory  measurements  program)  is 
calculated  in  this  report  to  be  minimal  for  the  beam  parameters  of  the  MICOM 
( S 3 )  laser  as  used  in  these  tests  and  therefore  is  not  expected  to  signifi¬ 
cantly  affect  the  calculations  of  the  clearing  phenomenon. 
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TABLE  /.  PARAMETERS  SHOWING  EFFICIENCY  OF  PUNCH-THROUGH  FOR  ETQTAL  =  40  J  FOR  CENTER  LOBE. 


On-Target 

Energy 

Without 


Aerosol  t'  Without 

Absorption  Time  to  Clearing, 

Coefficient  «  x  L  Evaporate  E/40 
(km-1)  (L  -  10-1km)  Aerosol  for  a  x  L 


t'  x  40 
~^pulse 


1.52  x  10-7 
1.54  x  10-7 
1.60  x  10-7 
1.68  x  10-7 


(5)  x  40 
(J) 


1.78 

X 

10-7 

'  0.86 

0.944 

1.91 

37.80 

1.87 

X 

10-7 

0.82 

0.942 

1.92 

37.70 

1.97 

X 

10'7 

0.78 

0.938 

1.92 

37.50 

2.08 

X 

10-7 

0.74 

0.933 

1.92 

37.40 

2.30 

X 

10-7 

0.6/ 

0.930 

1.93 

37.20 

2.56 

X 

1U"7 

0.61 

0.920 

1.95 

36.80 

4.34 

X 

10-7 

0.37 

0.864 

2.01 

34.60 

7.56 

X 

10"7 

0.22 

0.764 

2.08 

30.60 

1.41 

X 

10-6 

0.14 

0.559 

2.47 

22.40 

1.63 

X 

10"b 

0,12 

0.491 

2.45 

19.60 

1.92 

X 

10-b 

0.11 

0.400 

2.64 

16.00 

2.34 

X 

10-6 

0.10 

0.269 

2.93 

10.80 

2.70 

X 

10-b 

0.095 

0.156 

3.21 

6.24 

2.97 

X 

10-b 

0.093 

0.072 

2.45 

2.88 

>3.2 

X 

10-6 

0.093 

— 

---- 

— 

>3.2 

X 

10"b 

0.091 

— 

_ 

On-Target 
Energy 
With 
Clearing 
(6)  +  ( 7 ) / 40 
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APPENDIX  A 


CONSTANTS  USED  IN  THE  CALCULATIONS  OF  NONLINEAR  EFFECTS 

1.75  g/cc  (bulk  density) 

0.475  (cal/gm)/  C  (specific  heat  capacity) 

266  cal/gm  (latent  heat  of  vaporization) 

323  K  (ambient  temperature) 

443.7  K  (boiling  point  temperature) 


0.91  cm  center  spot 
0.778  cm  side  lobes 


(beam  radius) 


APPENDIX  C 


NEAR,  NEAR  FIELD  BURN  PATTERN 


7715  0.825  39.5  kV 

/ 7 16  0.900  36.0  kV 

7717  1.500  35.9 

7718  1.560  35.3 

7719  1.512  35.9 

7720  1.220 

1.320 
1.300 


7722 
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35.4 

35.4 

29.9  (prefire) 


APPENDIX  D 


APPENDIX  E 


METHOD  OF  ESTIMATING  ENERGY  DISTRIBUTION  IN 
THE  BEAM  FOR  LASER  BEAM 

Data  on  hand  consist  of: 

1.  Burn  pattern  reproductions  In  actual  size  for  field; 

Near  field  burn  pattern  (not  useful); 

3.  Some  dimensions  on  laser, 

that  is,  a.  Unstable  resonator  central  reflector  7.57  cm  "square"; 

b.  no  outside  limit  (beam  overfills  square  so  beam  energy 
cuts  off  at  about  10  cm.  This  is  not  adequate.  This 
outside  limit  must  be  calculated); 

c.  beam  take  off  and  reflector  supports  are  »  0.3 
inches  wide; 

d.  distance  from  laser  to  target  is  136  m. 

Thus,  the  output  aperture  geometry  is  as  shown  in  figure  E-l. 


Output  beam 
(intentional 
overfill ) 


Figure  E-l.  S3  laser  output  geometry. 


Approximation:  Use  central  arrays  of  orthogonal  twin  slit  patterns,  that  is, 


I 


sin2e 

e2 


(cos2  a]  , 


0  =  Y  b  sin  0  , 

n  J  „ 

a  -yd  sin  0  , 


where  b  is  the  effective  slit  width  at  the  output  of  the  unstable  resonator 
and  d  is  the  center-to-center  distance  of  the  output  beam  between  the  parallel 
sides. 

First  we  calculate  the  viewing  angle,  0,  in  terms  of  the  displacement  From  the 
center  of  the  arrays. 


1.36  x  104 


We  use  the  far  field  burn  pattern*  and  the  double-slit  part  of  the  calculation 
to  find  the  critical  dimension,  b. 


a  =  y  d  sin  0  , 

where  d  =  7.56  cm  +  1.0  x  (gap  width  in  centimeters). 


The  cos2  (yd  - - - I  factor  gives  the  fine  grain  pattern  under  the 

\ A  1.36  x  104  / 

single  aperture  intensity  envelope.  Maxima  in  the  pattern  are  determined  by 


d[cos ( c ’ dx ) 1 
<3x 


-c'  sin  c'dx 


0  , 


* 
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where 


- - - =  0.2179  cm'1 

1.36  x  104 


This  condition  is  met  when 


c'dx  =  0,  n,  2n  .  .  .  . 

First,  when  x  =  0  ; 

second,  when  x  =  — —  . 

c'd 

Then  from  the  photo  of  the  burn  pattern  (displacements,  x,  to  five  intensity 
lobes). 


d  =  <9.09  cm>  ±  0.07  cm  (or  ±8  percent)  . 

Calculating  the  quite  critical  parameter  b  for  the  diffraction  envelope 

sin2& 

~T~  » 

that  is,  9.09  cm  (twin  slit  center-to-center  separation) 

7,56  cm  (twin  slit  inside  dimension) 
b  =  1.53  cm  (2  x  (1/2)  x  slit  width)  . 

Therefore, 


j  b  =  4.535  x  103  ; 

£  d  *  2.694  x  104  . 

The  calculations  also  need  to  be  done  for  the  maxima  of  the  lobes.  The 
average  separation  from  the  determination  of  d  could  be  used, 

x  *  0,  1.6,  3.13,  4.77,  6.37  cm 

1.6  1.53  T64  1.6 
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that  is. 


x  =  <1.593>  ±  0.040  (or  ±2.5  percent) 


From  this, 


sin  e  =  1.171  x  10' 


Now  the  relative  intensities  in  the  lobes  of  the  far  field  pattern  can  be 
calculated  (table  E-l). 


TABLE  E-l.  CALCULATION  OF  RELATIVE  INTENSITIES  FOR  S3  LASER  BEAM  LOBES 


Can  See?  Spot  No.  Sin  e  Single  Slit  Double  Slit  Relative 


1.171 

X 

10‘4 

0.909 

2.342 

V 

10-4 

0.676 

3.513 

X 

10"4 

0.394 

4.684 

X 

10~4 

0.1604 

5.855 

X 

IQ"4 

0.0309 

7.026 

X 

10’4 

1.95  x  10-4 

o 

o 

2.0 

.999 

9.09  x  10"1 

,  999 

6.76  x  10'1 

.998 

3.93  x  10"1 

.997 

1.60  x  10_1 

.995 

3.08  x  10"2 

.994 

1.95  x  10-4 

It  will  be  assumed  that  all  the  energy  is  in  the  center  and  five  wing  spots  in 
each  direction  (21  spots). 


The  intensity  is  then  scaled  for  fraction  of  total  power  (10.675  norm). 


l? 

h 

O(ctr) 

1.87 

X 

10_1 

18.7 

1 

8.52 

X 

CM 

1 

O 

r-J 

8.5 

2 

6.33 

X 

IQ'2 

6.3 

3 

3.68 

X 

10"2 

3.7 

4 

1.50 

X 

10"2 

1.5 

5 

2.89 

X 

IQ-3 

0.3 

The  final  step  here  is  to  obtain  the  lobe  associated  numbers  giving  the 
fractional  power  densities.  To  obtain  these  numbers  divide  the  above  values 
by  1.90  for  all  but  the  center  lobe.  For  that  lobe,  divide  by  2.60  cm2.  The 
result  is  in  the  text. 
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APPENDIX  F 


HIGH-ENERGY  SINGLE-PULSE  PROPAGATION  ESTIMATES 


High-energy  laser  (HEL )  propagation  is  dependent  upon  a  variety  of  atmospheric 
conditions  and  beam  parameters.  Some  of  the  more  important  atmospheric 
factors  involved  are  absorption,  turbulence,  and  crosswinds.  Thermal  blooming 
and  gas  breakdown  thresholds  may  be  calculated  from  these  factors.1  In  this 
appendix  the  occurrence  of  these  two  effects  on  a  CO  HEL  beam  Is  formulated 
in  terms  of  limitations  cn  pulse  length.2  Assuming  a  fixed  total  pulse 
energy.  If  the  actual  pulse  length  of  the  laser  beam  Is  too  short,  then  the 
problem  of  gas  breakdown  becomes  Important  and  is  the  dominant  mechanism  that 
restricts  the  energy  fluence.  On  the  other  hand,  if  the  pulse  length  is  too 
long,  the  energy  fluence  on  target  drops  significantly  as  a  result  of 
transient  thermal  blooming.  These  two  effects  deleterious  to  HEL  propagation 
are  thus  minimal  for  a  range  of  pulse  lengths,  tp,  between  some  minimum  pulse 

length,  tm,  due  to  gas  breakdown  and  less  than  some  "saturated"  pulse  length, 

ts,  due  to  transient  thermal  blooming  effects;  that  is,  <  tp  <  t$. 


The  saturation  time,  defined  as  the  time  at  which  the  Instantaneous  peak 
Irradiance  at  the  target  has  dropped  to  about  10  percent  of  Its  Initial 
value,3  may  be  written  as  follows: 


where 


0.08 


i  Vocp 

I'VsSbs^p 


x2z 

DE  i/2 
P 


nQ  =  refractive  index 

p0  =  density 

Cp  =  specific  heat  at  constant  pressure 

nt  =  coefficient  of  index  change  with  respect  to  temperature 

cs  =  acoustic  wave  velocity 


1p.  G.  Gebhardt,  1976,  "High  power  laser  propagation,"  Appl  Opt,  15:1479 

2 S.  L.  Glicker,  1971,  "Propagation  of  a  10. 6pm  laser  through  a  cloud  including 
droplet  vaporization,"  Appl  Opt,  10:644 

3P.  B.  Ulrich,  1973,  "Requirements  for  experimental  verification  of  thermal- 
blooming  computer  results,"  J  Opt  Soc  Am,  63:897 
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aabs  =  absorption  coefficient 

f  *  correction  factor  that  accounts  for  finite  attenuation  and 
is  a  function  of  atZ,  the  optical  depth; 

that  is, 

f  =  — - - 

at  =  total  extinction  coefficient 

q  =  correction  factor  that  accounts  for  varying 
p  degrees  of  focusing; 

that  is, 

f  .  where  X  =  £ 

8X-  (X  -  1)*  Ad 

Aq  =  1/e  beam  radius 

A.  =  vacuum  1/e  beam  radius  at  the  range  Z  for  an  infinite 
Gaussian  source  including  diffraction  effects  and  poor 
beam  quality, 


a  z  -  1  +  e'atZ 


that  is 


A2 

i 


82  V 


+  (Ao 


- 

Tf) 


} 
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where  the  0  parameter  characterizes  the  beam  quality  of  the  source  in  terms  of 
its  far-field  or  focused  beam  radius. 

,  _  2n 

*  ~x" 

X  =  wavelength  of  the  source 
Z  =  range 

R  =  focal  range  of  the  beam 
D  =2  ST~  Aq  =  aperture  diameter 
Ep  -  total  pulse  energy 

The  minimum  pulse  length  limit,  t^  on  the  maximum  peak  target  irradiance  as 
imposed  by  the  gas  breakdown  threshold  is  simplified  as  follows: 


t 


m 


1.6 


02  ep 
»2Z2iBD 


e’“tZ 


9 


where  IBD  is  the  gas  breakdown  threshold  for  the  atmosphere  depending  upon 
many  factors  such  as  dust  aerosol  density,  laser  wavelength,  pulse  length,  and 
focal  snot  size  as  discussed  by  Morgan.4  For  our  calculations  we  assumed 
typical  values  for  the  atmosphere  of  n0  =  1,  p0  =  1  kg/m3,  cp  -  103  J/kg  -  K, 
-nt  =  10*6K"1,  c$  =  340  m/s,  and  IBB  =  107W/cm2,  which  is  appropriate  for  X  = 
10. Sum. 

The  following  table  is  a  compilation  of  ts  values  for  various  absorption 

coefficients  as  the  density  of  the  WP  cloud  decreases  with  time  after  initial 
dispersion. 


4C.  G.  Morgan,  1975,  "Laser-induced  breakdown  of  gases,"  Rep  Prog  Phys,  38:621 
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'MSSUELi 


*1 


Center  lobe  (0) 

Energy  =  40.0  J 

Effective  beam  diameter  aperture  =  2.57  cm 


aabs(km_1) 

f 

q 

t$(s) 

50 

0.321 

3.81 

3.33  x  10"6 

25 

0.506 

3.81 

3.75  x  10“6 

10 

0.736 

3.81 

4.92  x  10"6 

5 

0.852 

3.81 

6.46  x  10"6 

2.5 

0.922 

3.81 

8.79  x  10"b 

1.0 

0.967 

3.81 

13.6  x  10~6 

First  side  lobe  (1) 

Energy  =  18.1  J 

Effective  beam  diameter  aperture  =  2.20  cm 

aabs  1 ^ 

f 

q 

i  ( s ) 

50 

0.321 

7.36 

4. 

16 

lX> 

O 

r-* 

X 

25 

0.506 

7.56 

4. 

09 

X  10 '6 

10 

0.736 

7.36 

6. 

14 

10 

1 

o 

r -4 

X 

5 

0.852 

7.36 

8. 

08 

X  10“6 

2.5 

0.922 

7.36 

10 

.9 

U> 

1 

O 

X 

1.0 

0.967 

7.36 

16 

.9 

1 

o 

X 
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Second  side  lobe  (2) 

Energy  =  13.5  J 

Effective  beam  diameter  aperture  =  2.20  cm 


aabsUnf1) 

f 

q 

ts(s) 

50 

0.321 

7.36 

4.82  x  10“6 

25 

0.506 

7.36 

5.43  x  10~6 

10 

0.736 

7.36 

CD 

1 

o 

H 

X 

r — < 

• 

r-** 

5 

0.852 

7.36 

9.35  x  10‘6 

2.5 

0.922 

7.36 

12.7  x  10"6 

1.0 

0.967 

7.36 

19.6  x  10“6 

Third  side 
Energy  =  7 
Effective 

lobe  (3) 

.88  J 

beam  diameter  aperture  = 

2.20  cm 

aabs^m 

f 

q 

ts(  s ) 

50 

0.321 

7.36 

6.31  x  10'6 

25 

0.506 

7.36 

7.10  x  10-6 

10 

0.736 

7.36 

9.31  x  10"6 

5 

0.852 

7.36 

12.2  x  10-6 

2.5 

0.922 

7.36 

16.6  x  10~6 

1.0 

0.967 

7.36 

25.7  x  10"6 
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Fourth  side  lobe  (4) 

Energy  =  3.2  J 

Effective  beam  diameter  aperture  =  2.20  cm 


«abs(km"1) 

f 

q 

ts(s) 

50 

0.321 

7.36 

9.89  x  10"6 

25 

0.506 

7.36 

o 

H 

X 

0—4 

H 

r— < 

10 

0.736 

7.36 

14.6  x  10‘6 

5 

0.852 

7.36 

19.2  x  10~6 

2.5 

0.922 

7.36 

26.1  x  10"6 

1.0 

0.967 

7.36 

40.3  x  10‘6 

The  following  table  is  a  compilation 

absorption  coefficients  as  the  density 
time  after  initial  dispersion 

of  typical  tm  values  for  various 
of  the  WP  smoke  cloud  decreases  with 

D  (m)  Ep  (J) 

XBD 

(W/cm2 ) 

a 

(km-1 ) 

^m(s ) 

2.57  x  10“7  80 

107 

50 

5.07  x  10~9 

2.56  x  10“7  80 

107 

25 

6.17  x  10-8 

2.57  x  10'7  80 

107 

10 

2.77  x  lO"7 

2.57  x  10'7  80 

107 

5 

4.56  x  lO-7 

2.57  x  lO"7  80 

107 

2.5 

5.85  x  lO’7 

2.57  x  lO"7  80 

107 

1.0 

6.81  x  lO-7 
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Matching  the  average  pulse  length  tp  of  the  Sa  high  energy  CO^  laser  output  to 
the  minimum  pulse  length  requirement  and  to  the  upper  limit  on  the 
saturation  thermal  blooming  pulse  length  ts,  we  find  that  tp  Is  well  above  the 
gas  breakdown  threshold  for  a  variety  of  total  extinction  coefficients 
Also  for  absorption  coefficients  less  than  50  km"1,  tp  is  below  the  sate  ated 
pul  e  length  tg.  For  example,  at  the  center  lobe  (0)  with  high  aabs  anu  ~ 
50  km"1,  tp  ~  3.2ms  meets  the  requirements  to  maximize  the  single  pulse 
fluence  delivery  since  0.005ms  <  tp  <  3.33ms.  At  relatively  low  total 

extinction  values  ~  1.0  km"1,  tp  is  again  well  within  the  pulse  length 
limitations  imposed,  that  is,  0.58ms  <  tp  <  13.6ms  in  the  center  lobe  (0). 
Thus  we  see  that  thermal  blooming  and  gas  breakdown  effects  are  not 
significant  factors  in  punch-through  clearing  as  formulated  in  this  article 
for  the  S3  CO.,  HEL. 
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